Investigation of the relation between local diffusivity and local inherent structures in 

the Kob- Andersen Lennard-Jones model. 
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We analyze one thousand independent equilibrium trajectories of a system of 155 Lennard Jones 
particles to separate in a model-free approach the role of temperature and the role of the explored 
potential energy landscape basin depth in the particle dynamics. We show that the diffusion coef- 
ficient D can be estimated as a sum over over contributions of the sampled basins, establishing a 
connection between thermodynamics and dynamics in the potential energy landscape framework. 
We provide evidence that the observed non-linearity in the relation between local diffusion and basin 
depth is responsible for the peculiar dynamic behavior observed in supercooled states and provide 
an interpretation for the presence of dynamic heterogeneities. 

PACS numbers: 



Supercooled liquids are characterized by dynamics 
which take place on at least two well-separated time 
scales: a fast microscopic dynamics (associated to atomic 
or molecular motions at fixed liquid structure) and a slow 
structural relaxation time t q (which requires structural 
changes and diffusion processes) . If one considers the 
trajectory of the system in configuration space, such a 
separation of time scales can be visualized as the ex- 
ploration (on the microscopic time) of a finite region of 
the potential energy landscape (a basin), followed on a 
much slower pace by the exploration of distinct basins 
0- The potential energy landscape (PEL) formalism 
introduced by Stillinger and Weber Q provides a ther- 
modynamic description which builds on such a picture. 
The resulting expression for the liquid free energy re- 
quires information on the number of basins, their local 
minimum energy ers, and their shape. Analysis of nu- 
merical simulations, recently reviewed in Ref. has 
provided direct quantification of the statistical proper- 
ties of the landscape of several models of glass- forming 
liquids |3, 111 [D, |3, ISl [l0( • For the case of the Kob- Andersen 
Lennard-Jones model^, the number Q(ers)deis 
of distinct basins of energy depth between e/s and 
e/s + deis follows a Gaussian distribution p|, Q 
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where e is the total number of distinct basins for a 
system of N particles, E a is the energy scale of the 
distribution and a 2 is the variance; with <j, E a and 
a dependent only of the system number density. Of- 
ten, it is convenient to define the configurational entropy 
S con f(eis) = k>B log fl(eis), which for the Gaussian land- 
scape becomes a quadratic function of e/s 
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When fl (e/s) is known, the system partition function 
can be calculated^] (under the assumption of the e/s- 



independence of the basin anharmonicities) and predic- 
tions can be provided for the T-dependence of the ex- 
plored basin depth < eis{T) >, the probability P(e/s, T) 
of sampling a basin of depth e/s at temperature T and 
the configurational entropy S con f(T). 

The application of the Stillinger- Weber formalism to 
the analysis of numerical simulations provides an ap- 
pealing picture of the T-dependent exploration of the 
landscape and a convenient method to quantify the free 
energy of supercooled states. Several recent studies at- 
tempted, in different wa ys, to relate landscape properties 
to dynamics E Q [ll E HJ HI . The outcome 
of these studies strongly support the reasonable hope that 
in super-cooled states structural properties are strongly 
connected to dynamic ones and motivate the present de- 
tailed and essentially unbiased study of the connection 
between landscape properties and dynamics. 

In this Letter we study one-thousand independent runs 
of the well characterized Kob- Anderson model^J. We 
limit ourselves to the smallest system which can be sim- 
ulated with periodic boundary conditions without chang- 
ing the model's range of interaction, i.e. 155 particles at 
number density 1.2 12]. Due to the small system size, 
all quantities show large fluctuations, which we exploit 
to probe wide overlapping regions of e/s values under 
different T-conditions with the aim of disentangling the 
roles of e/s and T in the dynamics. 

We simulate six different T between 0.446 and 0.65. In 
this supercooled T- window, the decay of correlation func- 
tions clearly show the presence of two time scales and 
the diffusion coefficient D varies more than two orders 
of magnitude. For each T, we equilibrate each indepen- 
dent trajectory in the canonical ensemble and estimate 
P(e/s,7 1 ), which is found (see Fig. to be described 
rather well by a Gaussian distribution of T-independent 
variance a 1 = 18 ± 2, centered around a T-dependent 
mean < e/s > (T) ~ c 2 /T, in full agreement with the 
Gaussian landscape model and with previous results [2^] • 
The large a 2 guarantees a significant overlap of the land- 
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FIG. 1: Probability distribution P(eis,T) of exploring an IS 
of depth eis at temperature T. Black filled circles: equi- 
librium data for the N=155 system. Green full line: Gaus- 
sian distribution of variance 18. The figure also shows, with 
red filled squares the distribution P(< eis >i,T) (where 
< eis >i is not the instantaneous value but eis averaged 
over time £ms£>=i for the trajectory i and with a blue dashed 
line a Gaussian distributions with the same variance as P(< 
eis >i,T). 



scape region sampled at each T. 

Next we attempt to associate each of the independent 
configurations to an eis value characterizing the location 
of the system in the PEL and to a Di value characterizing 
its local diffusion coefficient. We further propagate each 
of the configurations for a time i<jwsD>=i, defined as 
the time at which the type-A particles mean square dis- 
placement (MSD), averaged over all trajectories and all 
A particles, reaches one. Different trajectories show very 
different dynamical properties. As an example, Fig. [2J-a 
represents the MSD, as a function of t for two differ- 
ent trajectories at the same T. Each trajectory i has a 
rather well defined diffusion coefficient Di which can be 
extracted from the slope of the MSD vs time. Fig. [SJb 
shows the i-evolution of eis f° r the same two trajecto- 
ries. The more mobile system is characterized by a larger 
and more fluctuating eis- Within t < MSD>=i, the set of 
eis values explored is far from covering the full range 
of values characteristic of that T. The system retains 
memory of the starting configurations and this behavior 
is enhanced at lower T. A quantification of the location 
of the ^-configuration on the PEL can thus be provided 
by < eis >»j defined as the value of eis averaged over 
the time interval from zero to t < MSD>=\- 

Fig. n also show the probability distribution P(< 
e is >i,T). While P(eis,T) are well approximated by 
a Gaussian distribution of average < eis{T) > and con- 
stant variance (in agreement with EqQ), P{< eis >*) are 
symmetric and well modelled by a Gaussian only at high 
T. At low T P{< eis >i,T) becomes strongly asymmet- 
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FIG. 2: a) MSD as a function of time for two different equili- 
brated trajectories at the same thermodynamical state point 
(T = 0.446, p = 1.2). The MSD is here averaged over all 
type-A particles of the system as well as over several starting 
times within the interval < t < t<MSD>=i- The dashed 
line represents the MSD averaged over all trajectories and it 
provides the definition of t<MSD>=i- b) Inherent structure 
energies as a function of time for the same trajectories shown 
in a). 

ric, and its variance increases on cooling. This originates 
in the different mobility of the system depending on the 
initial eis- Indeed, while systems with eis associated to 
the high energy tail of P(eis, T) manage to explore dif- 
ferent basins, systems in the low energy tail hardly lose 
memory of the initial configuration within £<msd>=i • In 
this respect, the observed asymmetry confirms that there 
is a strong coupling between < eis >i an d -C^, which we 
now try to quantify. 

Fig. 0-a shows Di as a function of < eis >i at several 
T. Data from different T overlap in a wide window of eis 
values. Data show that the local diffusivity is function of 
T not only via the T-dependence of eis but also directly, 
since for each chosen eis value, data for higher T always 
lay above data at lower T. The fact that basins of the 
same depth are sampled at different T allows us to evalu- 
ate the T dependence of D at fixed basin depth and sort 
out the relative contribution of T and eis to dynamics. 
To distinguish between the T and eis dependence of D 
we hypothesize that the law relating these two quantities 
has the form: 

D{ei Sl T) = fl^-^, (3) 

a form compatible both with the idea of activated dy- 
namics as well as with the Adam-Gibbs (AG) hypoth- 
esis ln(D/Doo) — T g B - (in which case l/f(eis) = 

Seonf( e Is)/ B), with B being a coefficient related to the 
minimum size of a cooperatively rearranging region. 

According to Eq. a plot of Tln(D{ei Sl T)/D 00 ) 
should produce a collapse of all data, independently of T, 
onto a master curve which can be identified as — /(ejg). 
Fig. 13-b confirms that indeed such a procedure gener- 
ates an impressive scaling of the data. We note that, in 
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FIG. 3: Diffusivity Di as a function of the corresponding 
average < eis >i- Different symbols indicate the six different 
simulated T. b) Tln(D(eis ,T) / Doo) as a function of eis- 
The solid line represent the fitting function ^ . 



this procedure, we use the value of — 0.155 obtain 
independently from high T simulations. If the result- 
ing f(eis) is interpreted according to activated models, 
one has to conclude that the height of the barriers in- 
creases significantly on cooling, diverging as the limiting 
eis value ei< = —7.81 is approached, a result which ap- 
pears to contradict, in the details, the analysis performed 
previously by Doliwa and Heuer jl9j . If f{eis) is inter- 
preted according to the AG equation, then l//(e/s) has 
to be proportional to S con f(eis), a quantity which has 
been previously calculated for this model using thermo- 
dynamic integration techniques 5, 23, 24j . 

The test between the two different estimates of S con f 
can be made even more stringent, if one realizes that, 
by fitting f(ejs) according to the relation /(e/g) = 
l IS - a2 )'i the fitting coefficient ag, a\ and a 2 can 
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Interestingly enough, since a 2 is known independently 
from the (T-independent) width of P(eis,T), ai pro- 
vides B, ao a and a 2 provides Eq (yielding respectively 
B = 1.85,a = 0.78,-Eo = -1140). Thus, an estimate of 
Sconf in absolute units is provided by the fit parameters. 
The comparison between Sconfobtained via thermody- 
namic integration (from Ref.|23j) and S CO nf obtained 
here from /(e/s) is shown in Fig. 0] The agreement be- 
tween the two independent estimates of S con f confirms 
the previously observed possibility of describing the dy- 
namics in this model according to the AG equation^] and 
to the interpretation of the landscape basins as statespH- 
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FIG. 4: Comparison between S CO nf evaluated from Eq. [2] 
and 21 (solid line), and S CO nf from I231 (dashed line). The 
symbols are the values obtained for each trajectory i. 



The combined use of the PEL thermodynamic formal- 
ism and of the relation between thermodynamic and dy- 
namics provided by Eq. |3| makes it possible to evaluate 
the T-dependence of D for the present model, — inde- 
pendently from the interpretation in terms of activated 
processes with eis dependent barrier or in term of AG 
— according to 



D(T) = J D(ei S ,T)P(eis,T)deis. 



(5) 



Fig. shows D evaluated directly from the MSD long 
time limit and from Eq. |S] The agreement between the 
two sets of independent data strongly support the pos- 
sibility of representing D as a sum over the local dif- 
fusivity of different parts of the landscape. This rep- 
resentation is very intriguing and offers an alternative 
way, within the landscape framework, to look at dy- 
namic heterogeneities, which have been widely studied 
theoretically, in simulations and experiments in recent 
years jH Hi |H HI HJ IH |H M, El IH E3? ? ]• 
Indeed, the nonlinear relation between D and eis (Eq.|3J) 
implies that at low T, in equilibrium, the system can be 
represented by coexisting regions in space with very dif- 
ferent local mobility. As seen in Fig. |31 at the lowest 
investigated T, the Di values cover about three differ- 
ent orders of magnitude. It is important to stress that 
the inhomogeneities in the dynamics are, in the present 
case, strongly associated to structural properties and in 
particular to the depth of the IS locally sampled. 

This extremely wide distribution of Di values suggests 
also a possible explanation of the T-dependence of the 
product Dr a upon supercooling, analyzed with reference 
to heterogeneous dynamics in, e. a . l2rj I27L l30t l3ll 
13^ . along the lines proposed in Ref. |26j|. Indeed, when 
a distribution is wide, the mean and the inverse of the 
average of its inverse can be rather different, since they 
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FIG. 5: a) Diffusivity constant as from data compared with 
the diffusivity constant evaluated from Eq. |K| b) Product 
< D X 1/D > form data as a function of T. The dashed 
line is a guide for the eye. 



weigh different parts of the distribution. To evaluate this 
effect in the present case we compare in Fig. the T 



dependence of the product < D >< 1/D >. On cooling, 
the product starts to deviate from one, consistently with 
the observations for the model liquid studied here that 
the product Drj/T |3(J and -DT a [37| starts to increase on 
cooling. 

In summary, by considering large numbers of realiza- 
tions of the trajectory of a small system of particles, we 
have shown that the dependence of the diffusion coeffi- 
cients for each such trajectory on the depth of the po- 
tential energy minima explored can be extracted in a 
model-free manner, and that the average diffusion co- 
efficient for the system can be estimated as a sum over 
contributions of the sampled basins, establishing a con- 
nection between thermodynamics and dynamics in the 
potential energy landscape framework. We provide evi- 
dence that the observed non-linearity in the relation be- 
tween local diffusion and basin depth is responsible for 
the peculiar dynamic behavior observed in supercooled 
states and provide an interpretation within the energy 
landscape picture for the presence of dynamic hetero- 
geneities. 
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thanks the University of Roma for hospitality. 
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